Nature of 3D Bose Gases near Resonance 



(N 

O 

(N 

in 



C/3 

I 

^ ■ 



O 

>: 

oo 

d 



X: 



Dmitry Borzov^, Mohammad S. Mashayekhi^, Shizhong Zhang^, Jun-Liang Song'^ and Fei Zhou^'^ 
Department of Physics and Astronomy, The University of British Columbia, Vancouver, B. C, Canada V6T1Z1^ 
Department of Physics, The Ohio-State University, Columbus, Ohio 43210, USA^ 
Institute for Quantum Optics and Quantum Information of the Austrian Academy of Sciences, 6020 Innsbruck, Austria^ 
Institute for Advanced Studies, Tsinghua University, Beijing 100084, Chinch 
(Dated: February 7, 2012) 

In this paper, we explore the nature of three-dimensional Bose gases at large positive scattering 
lengths via resummation of dominating processes involving a minimum number of virtual atoms. We 
focus on the energetics of the nearly fermionized Bose gases beyond the usual dilute limit. We also 
find that an onset instability sets in at a critical scattering length, beyond which the near-resonance 
Bose gases become strongly coupled to molecules and lose the metastability. Near the point of 
instability, the chemical potential reaches a maximum, and the effect of the three-body forces can 
be estimated to be around a few percent. 



I. INTRODUCTION 

Recently, impressive experimental attempts have been 
made to explore the properties of Bose gases near Fesh- 
bach resonance [H-Q- In these experiments, it has been 
suggested that when approaching resonance from the side 
of small positive scattering lengths in the upper branch, 
Bose atoms appear to be thermalized within a reason- 
ably short time, well before the recombination processes 
set in, and so form a quasistatic condensate. Further- 
more, the life time due to the recombination processes is 
much longer than the many-body time scale set by the 
degeneracy temperature. This property of Bose gases 
near resonance and the recent measurement of the chem- 
ical potentials for a long-lived condensate by Navon et 
al. [3| motivate us to make further theoretical investi- 
gations on the fundamental properties of Bose gases at 
large scattering lengths. 

The theory of dilute Bose gases has a long history, 
starting with the Bogoliubov theory of weakly interact- 
ing Bose gases [1] . A properly regularized theory of dilute 
gases of bosons with contact interactions was first put 
forward by Lee, Huang, and Yang [5] and later by Be- 
liaev 01 , who developed a field-theoretical approach. 
Higher-order corrections were further examined in later 
years Q . Since these results were obtained by apply- 
ing an expansion in terms of the small parameter V na^ 
(here n is the density and a is the scattering length), it is 
not surprising that, formally speaking, each of the terms 
appearing in the dilute-gas theory diverges when the scat- 
tering lengths are extrapolated to infinity. As far as we 
know, resummation of these contributions, even in an ap- 
proximate way, has been lacking [lo| . This aspect, to a 
large extent, is the main reason why a qualitative under- 
standing of Bose gases near resonance has been missing 
for so long. 

There have been a few theoretical efforts to under- 
stand the Bose gases at large positive scattering lengths. 
The numerical efforts have been focused on the energy 
minimum in truncated Hilbert spaces, which have been 
argued to be relevant to Bose gases studied in experi- 



ments prl - fisj . These efforts are consistent in pointing 
out that the Bose gases are nearly fermionized near res- 
onance. However, there are two important unanswered 
questions in the previous studies. One is whether the en- 
ergy minimum found in a restricted subspace is indeed 
metastable in the whole Hilbert space. The other equally 
important issue is what the role of three-body Efimov 
physics in the Bose gases near resonance is. 

Below we outline a nonperturbative approach to the 
long-lived condensates near resonance. We have applied 
this approach to explore the nature of Bose gases near 
resonance and to address the above issues. One concept 
emerging from this study is that a quantum gas (either 
fermionic or bosonic) at a positive scattering length does 
not always appear to be equivalent to a gas of effectively 
repulsive atoms; this idea, which we believe has been 
overlooked in many recent studies, plays a critical role in 
our analysis of Bose gases near resonance. 

Our main conclusions are fourfold: (a) energetically, 
the Bose gases close to unitarity are nearly fermionized, 
i.e., the chemical potentials of the Bose gases approach 
that of the Fermi energy of a Fermi gas with equal mass 
and density; (b) an onset instability sets in at a positive 
critical scattering length, beyond which the Bose gases 
appear to lose the metastability as a consequence of the 
sign change of effective interactions at large scattering 
lengths; (c) because of a strong coupling with molecules 
near resonance, the chemical potential reaches a maxi- 
mum in the vicinity of the instability point; (d) at the 
point of instability, we estimate, via summation of loop 
diagrams, the effect of three-body forces to be around a 
few percent. 

Feature (a) is consistent with previous numerical calcu- 
lations |llHl3| : both (b) and (c) are surprising features, 
not anticipated in the previous numerical calculations or 
in the standard dilute- gas theory [B','?']. Our attempt here 
is mainly intended to reach an in-depth understanding of 
the energetics, metastability of Bose gases beyond the 
usual dilute limit as well as the contributions of three- 
body effects. The approach also reproduces quantitative 
features of the dilute-gas theory. In Sec. H and Ap- 
pendixes A-C, we outline our main calculations and ar- 
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guments. In Sec. Ill, we present the conclusion of our 
studies. 



II. CHEMICAL POTENTIAL, METASTABILITY 
AND EFIMOV EFFECTS 



The Hamiltonian we apply to study this problem is 
H = ^(ek - /i)&lc&k + 2Uono ^ 6[6k 

k k 
k k 

+ ;^V^Xl^q^k'+a6_k'+a + H.c. 

k'.q 



2n 



II ^I+a&-k+fW§^-k'+H +H.C. (1) 

k,k',q 

Here ek = |k|^/2TO, and the sum is over nonzero mo- 
mentum states. Uq is the strength of the contact in- 
teraction related to the scattering length a via Uq^ — 
m(47ra)^^ — f^^^ X]k(2^k)~^, and is the volume, no 
is the number density of the condensed atoms and fi is 
the chemical potential, both of which are functions of a 
and are to be determined self-consistently. The chemical 
potential fi can be expressed in terms of i?(no,/i), the 
energy density for the Hamiltonian in Eq. (1), with hq 
fixed 



dE{nQ,fj.) 
dno 



,E{no,^i) = gM{no,^i)^, (2) 



M=2 



where gniM = 2,3, ...) are the irreducible M-body po- 
tentials that we will focus on below. The density of con- 
densed atoms no is further constrained by the total num- 
ber density n as 



no 



dE{nQ, ^) 
9/i 



(3) 



In the dilute limit, the Hartree-Fock energy density is 
given by Eq. ©, with g2 = Aira/m and the rest of the 
potentials gM,M = 3,4... set to zero. The one-loop con- 
tributions to gM for M = 3,4, ... in Figs. 1(c) and 1(d) 
all scale like 32 Vna-^, and their sum yields the well-known 
Lee-Huang- Yang (LHY) correction to the energy den- 
sity @. When evaluated in the usual dilute-gas expan- 
sion, g2 as well as one-loop contributions formally diverge 
as a becomes infinite. Below we regroup these contribu- 
tions into effective potentials 32,3... at a finite density hq 
via resummation of a set of diagrams in the perturbation 
theory. The approximation produces a convergent result 
for fj,. 

Before proceeding further, we make the following gen- 
eral remark. In the standard diagrammatic approach 
[l^ |. the chemical potentials can have contributions from 



diagrams with L internal lines, S interaction vertices, 
and X incoming or outgoing zero momentum lines, and 
X = 2S — L. For the normal self-energy (Sn) and 
the anomalous counterpart (S02) introduced by Beliaev, 
by classifying the diagrams Hugcnholtz and Pines had 
shown that, in general, the following identity holds [l4| in 
the limit of zero energy and momentum: /i = En — E02. 
Following a very similar calculation, we further find that 



S;ii("-o,m) = M + "o^, 
dno 



(4) 



where fi ~ dE{no, n)/dno. The equality in Eq. ^ is ef- 
fectively of a hydrodynamic origin. Following Eq. 
the speed of Bogoliubov phonons 4] Vs can be di- 
rectly related to an effective compressibility dno/dfi via 
mVg = Ell — = ''T'0-§^, where the first equality is 
due to the Hugenholtz-Pines theorem on the phonon 
spectrum |T5|. Note that hydrodynamic considera- 
tions had also been employed previously by Haldane 
to construct the Luttinger-liquid formulation for one- 
dimensional (ID) Bose fluids [1^. When na^ is small, 
Eq. (U) leads to the well-known result. En = 2/i. 

The self-consistent approach outlined below is mainly 
suggested by an observation that a subclass of one-loop 
diagrams [shown in Fig. 1(c)] yields almost all con- 
tributions in the LHY correction (see below and Ap- 
pendixes A and B). Resummation of these and their iV- 
loop counterparts can be conveniently carried out by in- 
troducing the renormalized or effective potentials (72,3 as 
shown in Figs. 1(a) and 1(b), where all internal lines 
represent, instead of the noninteracting Green's function 
GQ^{e,'k) = e — ek + /i + the interacting Hartree- 
Fock Green's function, G^^(e, k) — e — ey^ — En + fi + iS. 
This approximation captures the main contributions to 
the chemical potential in the dilute limit because the 
renormalization of two-body interactions is mainly due 
to virtual states with energies higher than fj, where the 
Hartree-Fock treatment turns out to be a good approxi- 
mation. The self-consistent equation for /i can be derived 
by estimating 52,3,... ('^•0; A^) diagrammatically (see exam- 
ples in Fig. 1). When neglecting 33,4,... potentials in 
Eq. ([2]), one obtains 



n = no + ^gl{no,n) 



d^k 9Eii/9no 
(27r)3 (ek-|-Eii-M)2' 
1 - 
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(5) 



Equations (U) and ([S]) can be solved self-consistently. 

We first benchmark our results with the LHY correc- 
tion or Beliaev's results for fj, by solving the equations 
in the limit of small na3. We find /i = ^noa(l -I- 
3-\/27r'\/noa3 + ...), and the number equation yields an 
estimate uq/ti — (1 — -^^V na3 -|- ...). The second terms 
in the parentheses are of the same nature as the LHY cor- 
rection. Comparing to Beliaev's perturbative result for 
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FIG. 1. (Color online) Diagrams showing contributions to 
the total energy E{no,fi). The dashed lines are for fc = con- 
densed atoms, thick solid internal lines in (a) and (b) are for 
interacting Green's functions G~^(e,k) = e — ek — Tjii + fi + iS, 
and thin solid lines in (c) and (d) are for noninteracting 
Green's function G(^^(e, k) = e — e^^ + fi + iS. (a) The blue cir- 
cle is for g2 {no , fJ.) ; vertices here represent the bare interaction 
Uo in Eq. (1). (b) (A'' = 1, 2, ...)-loop diagrams that lead to 
the integral equation for G3(— 377,p) in Eq. (7). Note that the 
usual tree-level diagram violates the momentum conservation 
and does not exist; the one-loop diagram has already been in- 
cluded in g2{no, /i) and therefore needs to be subtracted when 
calculating gs{no, /i). Arrowed dashed lines here as well as in 
(c) and (d) stand for outgoing condensed atoms, and the re- 
maining dashed lines stand for incoming ones, (c) and (d) The 
tree level and examples of one-loop diagrams that yield the 
usual Lee-Huang- Yang corrections in the limit of small na'^. 
The self-consistent approach contains contributions from (c)- 
type diagrams but not (d)-type ones (see further discussion 
in the text). Patterned green circles also represent the sum 
of diagrams in (a), but with thin internal lines, or the non- 
interacting Green's function Go lines. All vertices are time 
ordered from left to right. 



chemical potential, fi ~ ^noa{l + -^=y'noa^ + ...) @, 

and for the condensation fraction no/n = 1 — ^^-s/ na^ + 
one finds that the self-consistent solution reproduces 
99.96%(= 97rV2/40) of the Beliaev's correction for the 
chemical potential, and 83.30%(= 37r\/2/16) of the de- 
pletion fraction in the dilute limit. Technically, one can 
further examine (?2('t-0:M) by expanding it in terms of a 
and Ell and then compare with the usual diagrams in the 
dilute gas theory One indeed finds that g2{no,fi) in 
Eq. ([5]) effectively includes all one-loop diagrams with 
X — 3,4,5,... incoming or outgoing zero-momentum 
lines that involve a single pair of virtually excited atoms 
[between any two consecutive scattering vertices; Fig. 
1(c)]. The one-loop diagrams with X = 4,5,... incom- 
ing or outgoing zero-momentum lines that involve mul- 



(a) (b) 




0.1 0.2 0.1 0.2 



FIG. 2. (Color online) (a) Chemical potential fi in units of the 
Fermi energy ep and (b) condensation fraction as a function of 
Beyond a critical value of 0.18 (shown as circles), the 
solutions become complex, and only the real part of fi is plot- 
ted; the imaginary part of fi scales like ep^a/acr — 1)^^^ near 
Scr- (However, the sharp transition would be smeared out 
if the small imaginary part of G3 is included.) Dashed lines 
are the result of the Lee-Huang- Yang theory, thin solid blue 
lines are the solution without three-body effects (i.e. — 0). 
Thick solid red lines are the solution with gs included; the 
momentum cutoff is A = lOOn^''^. The inset is the relative 
weight of three-body effects in the chemical potential as a 
function of An~^^^ at the critical point. 



tiple pairs of virtual atoms [Fig. 1(d)] have been left 
out, but they only count for less than 0.04% of Beliaev's 
result [TtI. 

Following the same line of thought, one can also ver- 
ify that g2{no,iJ,) further contains {N — 2, 3, 4, ..)-loop 
contributions that only involve one pair of virtual atoms; 
each two adjacent loops only share one interaction vertex 
and are reducible. 173 (no, /i) included below, on the other 
hand, includes {N = 2, 3, 4, ..)-loop contributions with 
S — 4,5... interaction vertices that only involve three 
virtual atoms; two adjacent loops share one internal line 
instead of a single vertex [see Fig. 1(b)] and are irre- 
ducible, i.e., cannot be expressed as a simple product of 
individual loops. Effectively, we take into account all the 
virtual processes involving either two or three dressed ex- 
cited atoms in the calculation of the chemical potential 
fi by including the effective 52.3 (defined in Fig. 1) in 
Eq. (21). The processes involving four or more excited 
atoms only appear in (74.5... and are not included here; at 
the one- loop level following the above calculations, the 
corresponding contributions from the processes involving 
multiple pairs of virtual atoms are indeed negligible. 

A solution to Eq. ([S]) is shown in Fig. 2. An interesting 
feature of Eq. ([5]) is that it no longer has a real solution 
once Ti^/'^a exceeds the critical value of 0.18, implying an 
onset instability; this is not anticipated in the dilute-gas 
theory Q . This can also be illustrated by considering the 
two-body effective coupling constant 6*2 (Aq) as a func- 
tion of Ao [isj . a characteristic momentum that defines 



a low-energy subspace, 
G2(Ao) 



4tt 



(6) 



For Bose gases, it is appropriate to identify the relevant 
Aq as y^2mfi — A^. For positive scattering lengths, a 
not only defines the strength of interaction in the small 



A^ or dilute limit but also sets a 



scale for A^, above 



which the effective interaction becomes negative, i.e., 
G'2(A^) < if A^ > 7r/(2a). So as a approaches infinity, 
condensed atoms with a chemical potential fj, typically see 
each other as attractive rather than repulsive, resulting 
in molecules fl^. Thus, beyond the critical point the up- 
per branch atomic gases become strongly coupled to the 
molecules with a strength proportional to the imaginary 
part of fi. Consequently, we anticipate that /i decreases 
quickly beyond the critical scattering length due to the 
formation of molecules, leading to a maximum in /i in 
the vicinity of the critical point [l^ . 

A renormalization group approach based on atom- 
molecule fields was also applied in a previous study to un- 
derstand Bose gases near resonance [111, . Our results 
differ from theirs in two aspects. First, in our approach, 
an onset instability sets in near resonance even when the 
scattering length is positive, a key feature that is absent 
in that previous study. Second, when extrapolated to the 
limit of small na^, the results in Ref. [2T'| imply a cor- 
rection of the order of V na^ to the usual Hartree-Fock 
chemical potential but with a negative sign, opposite to 
the sign of LHY corrections. In a recent study jlSj, a 
self-consistent mean-field equation was employed, lead- 
ing to a similar conclusion as the approach in Ref. [2l| : 
the approach does not yield the correct sign of the LHY 
corrections. And so the onset instability pointed out in 
this paper, which is surprising from the point of view of 
dilute-gas theory, is also absent there. 

The chemical potential near the critical point can be 
estimated using Eq. ([S]) and is close to 0.9ei?, where 
ep = (67r^)^/'^?i^/'^/2m is the Fermi energy defined for 
a gas of density n. This is consistent with the picture of 
nearly fermionized Bose gases suggested by the previous 
calculations and experiments [s, 11 13, 21J. 

We now turn to the effect of 53 (ng, fi) on the chemical 
potential by including it in Eq. We estimate 173 by 
summing up all A^-loop diagrams with X — 3 incoming 
or outgoing zero momentum lines, which are represented 
in Fig. 1. All diagrams have three incoming or outgoing 
zero momentum lines but with N ~ 2,3,.. loops. The 
effect of three-body forces due to Efimov states [2^ was 
previously studied in the dilute limit [9]. The deviation 
of the energy density from the usual universal structures 
(i.e., only depends on na^ ) was obtained by studying the 
Efimov forces in the zero-density limit. The contribution 
obtained there scales like a^, apart from a log-periodic 
modulation [23|, and again formally diverges as other 
terms when approaching a resonance. There was also an 
interesting proposal of a liquid-droplet phase at negative 
scattering lengths but in the vicinity of a trimer-atom 
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FIG. 3. (Color online) ReG3(— 3?;, 0)»7 as a function of = 
Sii(no) — /i. Ea = A^/2m and A is the momentum cut- 
off. The imaginary part of G3 (not shown) is zero once 877 > 
l/(ma^). 



threshold (25|. 

It is necessary to regularize the usual a"* behavior at 
resonance in the three-body forces by further taking into 
account the interacting Green's function when calculat- 
ing the A^-Loop six-point correlators. Including the self- 
energy in the calculation, we remove the dependence 
that usually appears in the Bedaque- Hammer- Van Kolck 
theory for the three-body forces [2J|; when setting fi, Sn 
to zero, the equation collapses into the corresponding 
equation for three Bose atoms in vacuum, which was pre- 
viously employed to obtain the /3 function for the renor- 
malization flow in an atom-dimer field-theory model. 
The sum of loop diagrams in Fig. 1(b), G3{—3ri,p), sat- 
isfies a simple integral equation (m set to be unity; see 
Appendix C): 



G3(-3?/,p) = ^ J dqK{-3r];p,q) 



377 



3^' 



K[~3ri;p,q) = —In — -, 

pq p'^ + q'^ — pq + drj 



(7) 



where we have introduced rj = 'ShIuq) — fi. G3(— 377,0) 
is plotted numerically in Fig. 3. Three-body poten- 
tial 53(710, /i) is related to G3(— 377,0) via gsino,^) = 
6(7|ReG3(— 37/, 0) where G3 is obtained by further sub- 
tracting from G3 the one- loop diagram in Fig. 1(b) 
because its contribution has already been included in 
92(^0, p.). The structure of G3(— 37/,0) is particularly 
simple at a = -|-oo, as shown in Fig. 3: It has a de- 
sired log-periodic behavior reflecting the underlying Efl- 
mov states [2^ . When 377 is close to an Efimov eigenvalue 
B„ ~ -Bo exp(— 27r7i/so) [n = 1, 2, 3....,exp(27r/so) = 515] 
that corresponds to a divergence point in Fig. 3, the 
three-body forces are the most significant. When 37/ is 
in the close vicinity of zeros in Fig. 3, the three-body 
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forces are the negligible and Bose gases near resonance 
are dictated by the 172 potential. 



When including the real part of (73 (no, /x) in the calcula- 
tion of E(nQ, fj,), we further get an estimate of three-body 
contributions to the energy density and chemical poten- 
tial fi. The contribution is nonuniversal and depends on 
the momentum cutoff in the problem. For typical cold 
Bose gases, it is reasonable to assume the momentum 
cutoff A in the integral equation Eq. (O to be lOOn^/'^ or 
even larger. Quantitative effects on the chemical poten- 
tial are presented in Fig. 2. 



III. CONCLUSIONS 

In conclusion, we have investigated the energetics of 
Bose gases near resonance beyond the Lee-Huang- Yang 
dilute limit via a simple resummation scheme. We have 
also pointed out an onset instability and estimated three- 
body Efimov effects that had been left out in recent the- 
oretical studies of Bose gases near resonance 
Within our approach, we find that the three-body forces 
contribute around a few percent to the chemical potential 
and that the Bose gases are nearly fermionized before an 
onset instability sets in near resonance. 
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Appendix A: Solving Self-consistent Equation (5) in The Dilute Limit 



We apply Eq. (5) to calculate the leading-order correction beyond the mean-field theory. We notice that the 
equations for (72 and fi are arranged in a way that the next-order correction can be obtained by applying the results 
from the lowest-order approximation to the right-hand side. In the lowest-order approximation, we find Sn = 87rnoa 
and fi = 47moa; this leads to a correction to 52 as 



52 — 47ra 



{Ana) 



(Pk 



1 



1 



(27r) \<^k Cfc + Ai 



= Ana 



Similarly, from the relation 



driQ 



and ^^JJ- — 0, we can get the correction for the chemical potential /i as, 



{Airafnl f d^k 



II = Anano 
and the depletion fraction 



(27r)' (efe + fif 



Airano (1 + 3V27rnoa3) 



{2nf (cfe + lif 



TT 



For a comparison we list the results from the dilute-gas theory. 



ATTTloa 



40 



(Al) 

(A2) 
(A3) 

(A4) 
(A5) 
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Our self-consistent approach produces (= 99.96%) of Beliaev's result for the chemical potential, and (= 83%) 
for the depletion fraction. 



Appendix B: A Comparison Between the Self-Consistent Approach and the Dilute- Gas Theory 

hi the following, we show explicitly that our self-consistent equation corresponds to a subgroup of diagrams [in Fig. 
1(c)] in the usual dilute gas theory. The two-body T-matrix used in the dilute-gas theory [represented by the green 
circles in Figs. 1(c) and 1(d)] are obtained using the non- interacting Green's function G~^(e, k) = e — ek + + iO+; 
in the dilute limit, we can expand the T-matrix as 



t{u>, Q) = Ana 



1 -I- Ana 



(2^ 



(Bl) 



where w and Q are the total energy and momentum of the incoming atoms. The contribution from the first two 
diagrams in Fig. 1(c) are 



E(ci) 



Inaur. 



1 -f Ana 



1 



(c2) 



i2(0,0)ng f (fk 



(27r) 



1 



2nann 



(Ana) (167rnoa) 



dPk 



1 

(27r)^ \-k^ + 2^i + iQ+^ 

2 

2not{fi - efc,0) 



1 



{2Trf \-k^ + 2n + i0+ 



(B2) 
(B3) 



For the leading-order correction beyond the mean-field theory, it suffices to set i(/i — e^, 0) ~ Ana in Eq. (jB3l) and in 
higher-order diagrams. Similarly, we can get the contributions from the higher-order diagrams in this series, and the 
sum is 



2nan 



1 -I- Ana 



dPk 



{2nr 

2nan^{Ana) (167rano) 



1 



1 

fc2 + 2^ + iO+ l2 

d^k / 1 



{2nf V-fc^ + 2At-l-iO+ 



m-t-l 



2nanl{ 



1 -I- Ana 



(27r) V-fc^ + 2Ai-167rnoa k"^ 



(B4) 
(B5) 



We see that the energy given by the diagrams in Fig. 1(c) is exactly the same as the one used in our self-consistent 
equation, e.g., g2nf^/2, where g2 should be expanded as Eq. (|A1[) in the dilute limit. 

Next, we can sum up the rest of the one-loop diagrams that are not included in the self-consistent equations; they 
represent the lowest-order contributions to four- and six-body forces and so on. In the dilute limit, these diagrams 
[as shown in Fig. 1(d)] can be summed as 



E{d) = - (Anana) 



1 (2m -2)! 



Anarif) 



d^k 

{2nf 2m!(m- 1)! V 2efc - 2/i + 16no7ra 



2m-l 



(B6) 



Indeed, we can recover Beliaev's result by summing up one- loop diagrams in Figs. 1(c) and 1(d) as 



— (£'(c) + E(d)) = Annoa + Ana 



d^k 
J2nf 



(efe - ^ + QnriQa) 



Ann^a 



. 40/1 



\J (efc - /i -I- 87rano)^ - {AnariQ)' 



- 1 



Annaa 



(B7) 
(B8) 
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Appendix C: Including Three-body Forces In the Self-Consistent Equations 



We now calculate the amplitude of three-body scatterings corresponding to the processes described in Fig. 1(b). 
First, we consider a general case where three incoming momenta, instead of being zero, are ki = p/2 — q, k2 = p/2 + q, 
and = — p, and the outgoing ones are k'^ = p'/2 — q', kj = p'/2 + q', and kg = — p'. The scattering amplitude 
between theses states is then given by A{E — 877; p, p'), which represents the sum of diagrams identical to Fig. 1(b) 
except that the external lines carry finite momenta. 

For the estimate of three-body contributions of 173, we first treat the sum of diagrams in Fig. 1(b) as the limit 
of A{E — 877; p, p') when p and p' approach zero and the total frequency E is set to zero. It is therefore more 
convenient to work with the reduced amplitude G^{E — 877; p) = A{E — St]; p, 0), where p' is already taken to be zero. 
Gz{E — irj] p) itself obeys a simple integral equation, as can be seen by listing the terms in the summation explicitly. 
Indeed, when E is further set to zero, we find that the diagrams in Fig. 1(b) yield 

r ( \ 2 /■ jr(-877;p,g)g2 ~1 

G3(-3?7,p) = - / dq (CI) 

/,^,^,ir(-877;,,g)g^/^(-8,;g,g0g-^ 1 ^ ^^^^ 
!'/^ + 377-iy|^^-i'^" + 3^ 

where K{—Sri;p^ q) is the kernel defined in the main text. The sum of the above infinite series leads to the following 
integral equation of G3 : 

n ( \ 2 K{-3T^;p,q)q^ -1 ^2 K{-3v;p,q)q^ ^ ^ ^ 

G3{-3r],p) = - / dq = + - / dq = G'3(-877,g). (C8) 

^lq2 + 3,^-lj' + iV ^|g2 + 377-i 

Note that G'3(— 37/,0) defined above includes a diagram [the leftmost one in Fig. 1(b)] that has already been included 
in g2. To avoid overcounting, we subtract the first diagram in Fig. 1(b) from G3 as 



93 = 6.92 Re 



G3(-877,0)-^/..^(-^'^^«'«)«^ 



(C4) 



Alternatively, one can also carry out a direct summation of the diagrams in Fig. 1(b). It leads to a result that 
numerically differs very little from the estimation obtained above via an asymptotic extrapolation. For instance, a 
direct evaluation of those diagrams yields 

2 f ^^ Ki-2iT,0,q )q^ -1 

, , K{-2TT,0,q)q^ K{^3lT,q,q')q'' -1 , 
dqdq . ,2,0 + • • • (^6) 



G3 (-377, 0) = - / dq / " — (C5) 



The only difference between Eqs. (jC6p and (|C2I) is that the frequencies appearing in the first kernel K{E; 0, q) in the 
integrands and in the last denominators are now —277 instead of —877. 

One can easily verify that the sum can be written in the following compact form: 

i^(-277;0^)g2 ^ -1 
377 



G3(-37;,0) = - / dq ; " '^^^ -[^^+^3(^877,9)], (C7) 
_ I j_ 9^ + 277 



where G^{—Sri,p) is a solution of the following integral equation: 



2 r K{-37j;p,q)q 



2 



G3(-377,p) = - / dq / "-'■-^■- -[^^ + G3(-377,g)] (C8) 



Note that G^{—3ri,p) defined here describes an off-shell scattering between three incoming atoms with momenta 
p/2 — q, p/2 + q, and — p and three condensed atoms. And G3(— 37;,0) is not equal to G3(— 877,0), a consequence of 
the Hartree-Fock approximation we have employed here. G3(— 877,0) and G3 (—877,0) can be obtained numerically. 
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Finally, after subtracting the leftmost one- loop diagram in Fig. 1 (b) we again find the three-body contribution is 



93 



G3(-377,0)-- 

TT 



dq 



K{-2T^;0,q)q^ 



-1 



377 



1 g2 + 2?7 



(C9) 



Now we can include the three-body forces in a set of differential self-consistent equations similar to Eq. (5). We 
solve the equation numerically, and the results are shown in Fig. 2, where in the inset we show the momentum cutoff 
A dependence in the chemical potential. In our numerical program, we further use the approximation — 2g2, 

^ = 0, and Sii = (/3 = 2) to simplify the numerical calculations. We have tested other types of approximation 
schemes for the self-energy, such as Sn = Sirano or Sn — 2g2nQ. We find that the chemical potential and the value 
of the critical point na^^ are insensitive to approximation schemes. 
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